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This talk reviews the constraints imposed by binary-pulsar data on gravity theories, 
focusing on "tensor-scalar" ones which are the best motivated alternatives to general 
relativity. We recall that binary-pulsar tests are qualitatively different from solar-system 
experiments, because of nonperturbative strong-field effects which can occur in compact 
objects like neutron stars, and because one can observe the effect of gravitational radia- 
tion damping. Some theories which are strictly indistinguishable from general relativity 
in the solar system are ruled out by binary-pulsar observations. During the last months, 
several impressive new experimental data have been published. Today, the most con- 
straining binary pulsar is no longer the celebrated (Hulse- Taylor) PSR B1913-I-16, but 
the neutron star-white dwarf system PSR J1141— 6545. In particular, in a region of 
the "theory space" , solar-system tests were known to give the tightest constraints; PSR 
J1141— 6545 is now almost as powerful. We also comment on the possible scalar-field 
effects for the detection of gravitational waves with future interferometers. The presence 
of a scalar partner to the graviton might be detectable with the LISA space experi- 
ment, but we already know that it would have a negligible effect for LIGO and VIRGO, 
so that the general relativistic wave templates can be used securely for these ground 
interferometers. 



1. Introduction and solar-system constraints 

The most efBcient way to test a theory is to contrast its predictions with alternative 
models. Instead of just confirming or ruling out a particular theory, this method 
allows us to understand what features have actually been tested, and what kind of 
observations could be performed to test the remaining features. 

The best known example of such an embedding of general relativity in a space 
of alternatives is the so-called "parametrized post-Newtonian" (PPN) formalism^ 
which describes all possible metric theories of gravity in weak- field conditions, at 
order 1/c^ with respect to the Newtonian interaction. The basic idea was formulated 
by Eddingtonj^ who wrote the usual Schwarzschild metric in isotropic coordinates, 
but introduced some phenomenological parameters /3^^^ and 7^^^ in front of the 
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different powers of the dimensionless ratio Gm/rc?: 

9«='S«(l + 27-™|2)+o(l). (lb) 

General relativity corresponds to (3^™ — 7^^^ — 1, and is in perfect agreement 
with solar system experiments. At the time of the 10th Marcel Grossmann Meeting 
(MGX), the tightest published bounds on these parameters wer 

<2x 10-3, <6x 10-4. (2) 

An unpublishec3 stronger constraint was also known, [7^^^ — 1| < 4 x lO^"*, but 
an impressive new resultf^ has been released two months after MGX: 

7P™-l = (2.1±2.3)xlO-^ (3) 

Such bounds tell us that general relativity is basically the only theory consistent 
with experiment at the first post-Newtonian order. However, they do not constrain 
higher order terms in metric 1^, and the correct theory of gravity might differ 
significantly from general relativity in strong field conditions. Indeed, if R denotes 
the radius of a body, the ratio Gm/Rc^ « 10^^ for the Earth and « 10~® for the 
Sun, but it reaches ~ 0.2 for neutron stars, not far from the theoretical maximum 
of ^ for black holes. Pulsar observations can thus be used to test the strong-field 
regime of gravity. 

2. Binary-pulsar tests 

Experiment tells us that isolated pulsars are very stable clocks. A binary pulsar 
is thus a moving clock, the best tool that one could dream of to test a relativistic 
theory. Indeed, as the pulsar moves around its companion, the longitudinal Doppler 
effect modifies its observed spin rate, and we get a stroboscopic information on its 
orbital velocity. Several parameters characterizing its Keplerian orbit can then be 
extracted from an analysis of the pulse arrival times. For instance, the time between 
two maxima of the pulse frequency gives a measure of the orbital period P. One 
can also extract the eccentricity e, the longitude of periastron w, and the projected 
semimajor axis x along the line of sight {x = ^ smi, where i denotes the inclination 
of the orbit with respect to the plane of the sky). 

If high enough precision is achieved, relativistic corrections can be measured. In 
particular, if a system is observed long enough, the time derivatives of the Keple- 
rian parameters become available. In the case of the famous Hulse- Taylor binary 
pulsar PSR BIQIS+IG, three such "post-Keplerian" observables have been de- 
termined with great accuracy P (i) the Einstein time delay parameter jt, which 
combines the second-order Doppler effect (cx v\/2c'^, where va is the pulsar's ve- 
locity) together with the redshift due to the companion (cx Gms /rABC^ where tub 
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is the companion's mass and tab the pulsar-companion distance); (ii) the perias- 
tron advance u; (relativistic effect of order w^/c^); and (iii) the rate of change of 
the orbital period, P, caused by gravitational radiation damping (an effect of order 
/c' in general relativity, but generically of order /c^ in alternative theories, 
see below). In any theory of gravity, these three quantities can be computed in 
terms of the two unknown masses of the pulsar and its companion, and the equa- 
tions predictions(mA,mB) = observed values define three curves in the mass plane 
(or rather three strips if one takes into account experimental uncertainties). This 
provides 3 (observables) — 2 (unknown masses) — 1 test of the theory. If the three 
strips meet in a small region, there exists a pair of masses consistent with all three 
observables, and the theory passes the test. If they do not meet, the theory is ruled 
out. The upper-left panel of Fig. displays this mass plane in the case of general 
relativity, which is nicely consistent with experimental data within la error bars. 
This Figure also displays the mass plane for three other binary pulsars which have 
been timed accurately. 

For PSR B1534+12, five post-Keplerian observables have been measuredl^ 
7t, (jj and P have been described above, whereas r = Gitib/c^ and s = sini 
denote respectively the "range" (global factor) and the "shape" of the Shapiro time 
delay. General relativity passes again the 5 — 2 = 3 tests (at the l.Scr level for P, 
whose measurement is spoiled by a poorly known Doppler contribution due to the 
acceleration of the system towards the center of the Galaxy) . 

The timing measurements of jt, ^ and P for PSR J1141— 6545 are brand- 
new,El and have been presented at this parallel session PTl of the MGX meeting 
(see M. Bailes' contribution to the present proceedings). This binary pulsar was 
discovered in 1999, and has been quickly recognized as a "strange" system, which 
triggered several detailed studies. The pulsar is young 1.4 Myr) and thereby 
nonrecycled, as indicated by its rather slow pulse period (0.4 s). The orbital period 
is short (P = 4h 45min), therefore large relativistic effects are expected. The 
companion is a white dwarf at the 90% confidence level,'^ but the eccentricity of the 
orbit (e = 0.17) is surprisingly large of such an asymmetrical system. Indeed, almost 
all known neutron star-white dwarf binaries have a vanishingly small eccentricity. 
The explanation is that the neutron star was formed after the white dwarf in PSR 
11141—6545.'^ Initially, the progenitor's mass was too small to evolve into a neutron 
star, but it accreted enough matter from its companion, and it finally exploded as 
a type Ib/c supernova, giving a momentum kick to the newborn neutron star. 
Before precise timing measurements could be performed, the Keplerian parameters 
of this system could be determined by an analysis of the pulsar scintillation.EH This 
phenomenon, caused by diffraction in the interstellar medium, happens to have a 
characteristic timescale of minutes for PSR J1141— 6545, and this is quick enough 
to extract orbital information. Indeed, the periodic variation of this scintillation 
timescale can be interpreted as a consequence of the orbital motion of the pulsar. A 
fit of the observational data gave two possible orbital solutions, -'^^ only one of them 
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Figure 1. Mass plane (m^ = pulsar, rrig = companion) in general relativity for the four precisely 
timed binary pulsars. The various strips are consistent at the la level with the timing observables 
which label them. In the lower two plots, the shaded region corresponds to | sini| > 1 and is thus 
excluded. 

predicting a pulsar mass ~ l.Sm© consistent with all other known neutron star 
masses (1.35 ± 0.04)771©. [The fact that it is the lightest known is consistent with 
the formation scenario summarized above.] The recent timing measurementJ^ of 
this system nicely confirm (within general relativity) this orbital solution, as shown 
in the lower-left panel of Fig. ^ Indeed, the intersection of the oj and 7t curves 
gives for the masses 77i^/mQ = 1.30 ± 0.02 and 7713/7710 = 0.986 ± 0.02. Kepler's 
third law implies the following relation between the inclination angle i, the masses, 
and the observed quantities P and x: 

(777^ + 7775)2 " [tJ ~g~- 

Therefore, once 777^1 and 777^ have been determined thanks to two precise enough 
post-Keplerian observables, this "mass function" |0J provides the inclination angle 
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of the orbit. One finds i > 75° at the la level, showing that the orbit in nearly edge 
on. In such a situation, the mathematical bound^ | sini| < 1 can be used as an extra 
constraint in the mass plane (m^,ms). The excluded region is shaded in Fig. ^ 
and we do verify that the intersection of the strips is very close to the limit sini = 1. 
This extra constraint will be quite useful below to exclude some alternative theories 
of gravity. For the above masses and the observed Keplerian parameters of the 
orbit, general relativity predicts an orbital period derivative P — —3.8 x 10~^^, 
nicely consistent with the observed value P°^^ — (— 4 ± 1) x 10^^'^, as illustrated 
also in the lower-left panel of Fig. ^ The la relative errors are still rather large, 
but since they scale as they should reach ~ 1% by 2010. We will see anyway, 

in Sec. 4 below, that the present precision is already extremely constraining for 
alternative theories, and that PSR J1141— 6545 is the best system available to test 
gravity in the strong- field regime. 

The timing measurements for PSR J0737— 3039 were releasecP^^five months 
after the MGX meeting, but this system is so interesting that I must mention it 
here. Like PSRs B1913-K16 and B1534-I-12, this is a double neutron star system, 
but its orbital period is so short (P = 2h 27min) that huge relativistic effects are 
expected. [This short orbital period also implies that the system will merge in about 
85 Myr, much quicker than any other known neutron star binary; this suffices to 
increase by a factor 10 the estimated merger rate in our Galaxyf^ and thereby 
our chances to observe a strong gravitational wave signal in the LIGO/ VIRGO 
interferometers.] The periastron advance could be determined in only a few days 
of observation, and its value uj — 17°/yr is much greater than any of the above 
binary pulsars (the largest being 5.3°/yr for PSR J1141— 6545). A careful analysis 
of the data taken during several months also provided three other post-Keplerian 
parameters,^! corresponding to the Einstein (77-) and the Shapiro (r and s) time 
delays. The measure of s = sini shows that the orbit is again almost perfectly 
edge on: i = (87 ± 3)°. Therefore, the mathematical bound |sini| < 1 can again 
be used as an extra constraint in the mass plane. The shaded region is excluded 
in the lower-right panel of Fig. But the greatest novelty of this system is that 
the companion has also been observed as a pulsar.IElof course, this second pulsar 
cannot be recycled itself (since only the first formed could accrete matter from its 
companion), therefore its pulse period is slow (2.8 s) and its observation is not 
precise enough to provide extra post-Keplerian observables. However, its Keplerian 
orbital data do suffice to give an extra test. Indeed, if xa and xb denote the observed 
projected semi-major axis of both pulsars, the ratio xa/xb = ms/mA gives a very 
precise measure of the mass ratio, independently of the theory (and valid at any 



^Note that the Shapiro post-Keplerian parameter s = sinj has not yet been measured precisely 
enough to provide another test with PSR J1141— 6545 (it should be available within a few months). 
We are here using only the mathematical fact that it cannot be greater than 1. 
''Contrary to PSR J1141— 6545, the parameter s is measured for PSR J0737— 3039, and the cor- 
responding solid curve should not be considered as the mere boundary of the shaded region. Its 
width is larger than Icr error bars. 
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order in powers of 1/c, at least in general relativity and tensor-scalar theories). 
Moreover, if one exchanges the roles of bodies A and B, the mass function and 
the bound | sini| < 1 give a second excluded region in the mass plane. This is the 
reason why the white region has got a V shape in the lower-right panel of Fig. ^ 
and why the shaded (excluded) region covers almost all the plane. This will be quite 
important to constrain alternative theories of gravity in Sec. 4 below. The orbital 
period derivative P of this double-pulsar system has not yet been determined, but it 
should be available with reasonable precision within a few months. It will of course 
provide an extra strip in the mass plane, and thereby an extra test of relativistic 
gravity in the strong-field regime. Geodetic precession has already been observed in 
PSRs B1913+16 and B1534-I-12. This confirms the existence of such an effect, and 
it can be used to map the emission beam of the pulsar (which gives us information 
about its structure), but this does not provide an actual test of the theory. On 
the other hand, the geodetic precession period is predicted to be of about 70 years 
for both pulsars in PSR J0737— 3039, within general relativity. Several years of 
observation should thus provide an extra test of the theory, which would not have 
been possible with the other known binary pulsars presented above. Let us finally 
mention that because this double-pulsar system is almost perfectly edge on, one 
observes eclipses of the recycled pulsar by its companion, and modulations of the 
companion's pulses caused by the energy flux of the recycled pulsar. Therefore, 
this system will allow us to probe pulsar magnetospheres. Again, this cannot be 
considered as a test of the gravity theory, but as a powerful new observatory of the 
physics of neutron stars. 

Besides the above four precisely timed binary pulsars, other systems do provide 
extra tests, even if one does not measure enough post-Keplerian parameters to de- 
termine accurately the two masses as in Fig. Indeed, one can use a statistical 
argument on the pulsar's mass ttia ~ 1.35mQ and the a priori arbitrary inclination 
angle i to predict a probable value for the companion's mass rriB thanks to the 
Keplerian "mass function" Then, a single post-Keplerian observable may be 
compared to the prediction of a theory. For instance, in the neutron star-white 
dwarf binary PSR B0655+64, only an upper bound has been obtained on the 
orbital period derivative P. It is an order of magnitude larger than the predic- 
tion of general relativity, therefore this theory is obviously consistent again with 
experimental data. However, this upper bound suffices to rule out a wide range of 
alternative theories of gravity, which generically predict a value for P several orders 
of magnitude larger (see Sees. 3 and 4 below). Alternative theories also predict 
various effects which vanish identically in general relativity, caused by violations 
of the strong equivalence principle, of local Lorentz invariance or of conservation 
laws. The experimental upper bounds on such effects therefore just constitute "null 
tests" for general relativity, but have the capability of constraining other theories. 
We will mention such an example in Sec. 4 below. 
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3. Tensor-scalar theories of gravity 

As shown is the previous section, several tests of gravity are available in the strong- 
field regime, and general relativity passes all of them with flying colors. We now 
wish to embed Einstein's theory into a class of alternatives, in order to understand 
better which features have been tested, and to compare the probing power of the 
various tests. A generalization of the PPN formalism to all orders in 1/c" would 
need an infinite number of parameters. It is much more efficient to restrict our study 
to the most natural class of alternatives to general relativity, namely "tensor-scalar" 
theories, in which gravity is mediated by a tensor field (g^ij/) together with one or 
several scalar fields (cp). This class of models is privileged for many reasons, (i) They 
are mathematically consistent field theories, and do not involve any negative-energy 
mode nor any adynamical field, (ii) The existence of scalar partners to the graviton 
is predicted by all unified and extra-dimensional theories, notably superstrings. 
(iii) Scalar fields play a crucial role in modern cosmology, notably to explain the 
accelerated expansion phases of the universe (inflation, quintessence), (iv) Tensor- 
scalar models are the only consistent massless field theories able to satisfy the weak 
equivalence principle (universality of free fall of laboratory-size objects), (v) They 
are the only known theories satisfying "extended Lorentz invariance", i.e., such 
that the gravitational physics of subsystems, influenced by external masses, exhibit 
Lorentz invariance. (vi) They explain the key role played by /3^^^ and 7^^^ in the 
PPN formalism (the extra 8 parameters introduced by Nordtvedt and WilP vanish 
identically in tensor-scalar theories) . (vii) They are general enough to describe many 
different deviations fr om g eneral relativity, but simple enough for their predictions 
to be fully worked out.l^ 

In this paper, we will further restrict our study to theories which involve a single 
scalar field, and which satisfy exactly the weak equivalence principle. Like in general 
relativity, the action of matter is given by a functional S'm['0m,5/i!/] of some matter 
fields V'm (including gauge bosons) and one second-rank symmetric tensor g^j,. The 
difference with general relativity lies in the kinetic term of g^^. Instead of being a 
pure spin-2 field, it is here a mixing of spin-2 and spin-0 excitations. More precisely, 
it can be written as — exp[2a((/3)]g^,y, where a{ip) is a function of a scalar field 
ip, and g^i, is the Einstein (spin 2) metric. The action of the theory reads thus 

' d^xV^ (R - 2g^''d^^d,^) + l^m^e^^^'^^gj . (5) 
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[Our signature is — h++, R is the scalar curvature of g^i,, and g its determinant.] 

Our discussion will now be focused on the function a(ip), which characterizes 
the coupling of matter to the scalar field. Let us expand it around the background 
value of the scalar field, which can be chosen to vanish without loss of generality: 

a{ip) = aQip + ^/3o'P^ H (6) 

The slope ag measures the coupling strength of the linear interaction between matter 
and the scalar field, /3o is the quadratic coupling constant of matter to two scalar 
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lines, etc. [A diagrammatic representation is given to label the axes of Figs. 01 and |H 
below.] General relativity corresponds to a vanishing function a{(p) = 0, and Jordan- 
Fierz-Brans-Dicke theory to a linear function a{ip) = aoip, with aQ — 1/ (2cjbd + 3). 
As shown below, interesting strong-field effects occur when (3o ^ 0, i.e., when a{ip) 
has a nonvanishing curvature. 

At the first post-Newtonian order (i.e., when measuring effects of order f/c^ in 
weak- field conditions) , the predictions of tensor-scalar theories depend only on the 
first two parameters ao and (3a. The effective gravitational constant between two 
bodies and the Eddington PPN parameters read: 



- G(l + ai) , 
P™-f = -2a2/(l + a2) 

1 ao/?oQ!o 

2 {1 + a, 



PPN 



f = 



2\2 



(7a) 
(7b) 

(7c) 



0^ 



[The factor comes from the exchange of a scalar particle between two bodies, 
whereas ao/?oQ^o comes from a scalar exchange between three bodies.] The solar- 
system bounds ©-(jSl therefore impose that both and Q!ol/3o| must be small. 
This implies that the scalar field must be linearly weakly coupled to matter. On 
the other hand, the quadratic coupling strength /3o is not directly constrained if 
is small enough, and its sign can also be arbitrary. [Note that a{ip) is a coupling 
function and not a potential for the scalar field, therefore a negative /3o does not 
spoil the stability of the field theory] 

At higher post-Newtonian orders l/c", a simple diagrammatic argument shows 
that any deviation from general relativity involves at least two factors ao, and has 
the schematic form 



deviation from G.R. = ar. 



Gm f Gm \ ' 



(8) 



where Aq, Ai, ... are constants built from the coefficients ao, /3oj ... of expansion 
© . Since is experimentally known to be small, we thus expect the theory to be 
close to general relativity at any order. However, some nonperturbative effects may 
occur in strong- field conditions: If the compactness Gm/Rc^ of a body is greater 
than a critical value, the square brackets of Eq. (O can become large enough to 
compensate even a vanishingly small aQ. To illustrate this, let us consider a model 
for which ao vanishes strictly, i.e., which is perturbatively equivalent to general 
relativity: There is strictly no deviation from general relativity at any order in a 
perturbative expansion in powers of 1/c. A parabolic coupling function a((p) = 
\I3qlp^ suffices for our purpose. At the center of a static body, the scalar field 
takes a particular value (fc, and it decreases as 1/r outside. The energy of such a 
scalar field configuration involves two contributions, coming respectively from the 
kinetic term and from the matter-scalar coupling function in action 10). As a rough 
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estimate of its value, one can write 



Energy ; 



mc 



y^g/2 

Gm/Rc^ 



(9) 



When Po < 0, this is the sum of a parabola and a Gaussian, and if the compactness 
Gm/Rc^ is large enough, the function Energy((/Jc) has the shape of a Mexican hat, 
see Fig. 12 The value (pc = now corresponds to a local maximum of the energy. It 
is therefore energetically favorable for the star to create a nonvanishing scalar field 
ipc and thereby a nonvanishing "scalar charge" a'{ipc) — /?oVc- This phenomenon 
is analogous to the spontaneous magnetization of ferromagnets. 




Figure 2. Heuristic argument to explain the phenomenon of "spontaneous scalarization" . When 
/3o < and the compactness Gm/Rc^ of a body is large enough, it is energetically favorable to 
create a local scalar field different from the background value. The body becomes thus strongly 
coupled to the scalar field. 



This heuristic argument has been verified by explicit numerical calculations, 
taking into account the coupled differential equations of the metric and the scalar 
field, and usin g va rious realistic equations of state to describe nuclear matter inside 
a neutron star.l^ The correct definition of the linear coupling strength between a 
compact body A and the scalar field reads aA = d\n ■mA/d(po. It is plotted in Fig. 13 
for the particular model /3q = —6. One finds that there exists indeed a "spontaneous 
scalarization" above a critical mass (whose value decreases as —/3o grows). On the 
other hand, if /3o > 0, both the above heuristic argument and the actual numerical 
calculations show that \aA\ < \ao\- In that case, one finds that neutron stars are 
even less coupled to the scalar field than solar-system bodies. 

The scalar charge aA enters the predictions of the theory in the same way as ao 
in weak-field conditions. For instance, the effective gravitational constant between 
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neutron star 




> mA/niQ 



maximum maximum 
mass in GR mass 



Figure 3. Scalar charge aA versus baryonic mass ruA, for the model a{(fi) = —3ip'^ (i.e., /3o = —6). 
The solid line corresponds to the maximum value of qq allowed by the Eqs. 0, and the dashed 
line to CfQ = 0- The dotted lines correspond to unstable configurations of the star. 



two bodies A and B reads 

G% = G(l + aAas) , (10) 

instead of Eq. H7a|) . Similarly, the strong- field analogues of the Eddington parame- 
ters 7^^^ and (3^^^ involve products of aA, as, Pa = daA/dipo and Pb, instead 
of ao and (3o as in Eqs. (TJ. Since the scalar charge a a ~ 0.6 in the model of Fig.|31 
one thus expects deviations by ~ 35% from some general relativistic predictions. 
Moreover, the quadratic coupling strength Pa can take very large numerical values 
near the critical mass, like the magnetic susceptibility of ferromagnets. Therefore, 
even larger deviations from general relativity are found when the mass of a neutron 
star happens to be close to the critical one. An even more spectacular phenomenon 
occurs for the quantity as^ln Jyi/i9(/3o, involved in the post-Keplerian observable'^ 
7t- Indeed, when Pq < 0, this quantity blows up as ao 0: Paradoxically, a 
theory which is closer to general relativity in weak-field conditions predicts larger 
deviations in the strong-field regime! On the other hand, when Pq > 0, strong-field 
predictions are even closer to those of general relativity than in the solar system. 

The post-Keplerian observable which is the most affected by the presence of a 
scalar partner to the graviton is the orbital period derivative P. Indeed, the energy 
flux carried out by gravitational waves is of the form 



r Quadrupole , „ / 1 
Energy flux = i ^ ^ I 77 



r Monopolc 



Dipolc Quadrupole 



hclicity 2 

-(^ 



(11) 



hclicity 

The first curly brackets contain the prediction of general relativity. The second ones 
contain the extra contributions predicted in tensor-scalar theories. The powers of 



'^This expression enters ■yx because the inertia moment of the pulsar, Ia, is modified by the 
presence of a companion at a varying distance. 
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1/c give the orders of magnitude of the different terms. In particular, the monopolar 
and dipolar helicity-0 waves are generically expected to be much larger that the 
usual quadrupole of general relativity. However, the scalar monopole has the form 

Monopole _ G \ d{mAaA) ^ dijUBaB) ^ ^ f ^ ^2) 
c c \ dt dt \c^ / J ' 

and it reduces to order 0{l/e') if the stars A and B are at equilibrium, dt{mAOiA) = 
0, which is the case for all binary pulsars quoted in Sec. 2 above. However, this 
monopolar term would be huge in the case of a collapsing star, for instance. The 
dipole has the form 

Dipole G fGABmAirLBY, \2 , f ^ \ n Q^ 

and is usually much larger that a quadrupole of order l/c'. For instance, in a 
pulsar- white dwarf binary, the pulsar's scalar charge a a may be of order unity, like 
in Fig. 131 above, whereas the weakly self-gravitating white dwarf has a very small 
scalar charge ub ~ aoi constrained by Eqs. l(2Jl and (|7b|) . On the other hand, in 
a double-neutron star system, one expects tua ~ tub and therefore a a ~ cub, so 
that this dipolar contribution is considerably reduced (but may still be large with 
respect to the usual quadrupole of general relativity). Indeed, a dipole is a vector 
in space, and two strictly identical stars do not define a preferred orientation. 

It should be noted that even within Brans-Dicke theory, i.e., for a linear coupling 
function a{ip) = aQip, the dipolar contribution does not vanish identically. 
Indeed, in that case, one can show that ua = Q;o(l — 2sa), where sa ~ Gm/R(? 
is the compactness of the star, so that {a a — as)^ = 4ag(sA — ss)^. The dipolar 
contribution is thus proportional to the experimentally small factor = l/(2a;BD + 
3), but it does not vanish, and can still be larger than the usual quadrupole if 

> w^/c^ ^ 10~^. Note also that because sa « Gm/Rc^, this dipole ((T^ 
is formally of order 0{G^ /c^), so that one could not have obtained it by a naive 
calculation at linear order in G. It is however crucial to take it into account, because 
it is numerically of order 0{G/c') since Gm/Rc^ « 0.2 for a neutron star. 



4. Experimental constraints on tensor-scalar theories 

We saw in the previous section that solar-system experiments probe only the first 
two coefficients ao and /?o of the matter-scalar coupling function JHl . On the other 
hand, the strong-field predictions do depend on the full shape of this function. 
In order to compare easily the different constraints, we will restrict our study to 
a strictly parabolic function a{(p) = aof + ■^Po'p'^- Figure 21 displays the known 
constraint's at the time of the previous Marcel Grossmann meeting MG9, three 
years ago. Since the physics does not depends on the sign of ao, only a half plane 
(|ao|j/3o) is drawn. Each point on this figure represents a tensor-scalar theory. 
General relativity is at the origin ao = Po = 0, Brans-Dicke theory is on the 



12 



matter 

O 

locol 

t LIGO/VIRGO 




Figure 4. Constraints on generic tensor-scalar theories imposed by solar-system experiments, 
classic binary-pulsar tests, and future detections of inspiralling binaries with laser interferometers. 
The hatched region is allowed by all the tests. The doubly hatched one is also consistent with 
VLBI data.4 



vertical axis Po — 0, and the horizontal axis ao = corresponds to theories which 
are perturbatively equivalent to general relativity (see Sec. 3 above). For each 
theory (|q;o|,/3o), the mass planes of the various binary pulsars can be plotted like 
in Fig. ^ and if the different strips do not have a common intersection, the theory 
is ruled out. The allowed theories lie under and to the right of the various curves. 

The thin solid line labeled "solar system" corresponds to the bounds and 
the unpublished limit provided by very long baseline interferometrjS (VLBI) is 
materialized by a thin horizontal line. As discussed in the previous section, they 
impose a very small matter-scalar linear coupling strength ag, but do not constrain 
the quadratic coupling Po if cto is small enough. 

On the other hand, binary pulsars impose Po > —4.5, even for a vanishingly small 
OfQ. As discussed in the previous section, this constraint is due to the spontaneous 
scalarization of neutron stars, which occurs when — /3o is large enough. Equations (0) 
allow us to rewrite this inequality in terms of the Eddington parameters (3^^^ and 
7^^^^, which are both consistent with 1 in the solar system. One finds 

oPPN _ 1 
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The singular (0/0) nature of this ratio vividly expresses why such a conclusion could 
net be obtained in weak- field experiments, and underlines that binary-pulsar tests 
are qualitatively different. 

This bound is mainly due to the Hulse- Taylor binary pulsar PSR B1913+16. 
For PSR B1534+12, whose orbital period derivative P is not very well known, we 
have plotted the constraints imposed only by the four post-Keplerian observables 
u>, 7t, r and s. This gives much weaker constraints than PSR Bf9f3-|-16, but this 
is nevertheless a very important test because it does not depend on the radiative 
structure of the theory. Taking into account the measured value of P removes the 
horn-shaped region at the top-left of the dashed line. 

The second most constraining system is the neutron-star white dwarf PSR 
B0655+64, although it has not been timed precisely enough to determine the two 
component masses, and that only an upper experimental bound on its orbital period 
derivative P is known. Indeed, even by choosing conservative values for the two 
masses, such an asymmetrical system loses too much energy through dipolar scalar 
waves in most scalar-tensor theories, see Eq. (I13|l . 

Equation (|10|l above tells us that the acceleration of a neutron star A towards 
the center C of the Galaxy is proportional to (1 -I- aAOic), whereas a white dwarf 
B is accelerated with a force oc (1 + a sac)- Since a a as in general, there is a 
violation of the strong equivalence principle (SEP). This causes a polarization of 
the orbit of a neutron star- white dwarf system towards the Galaxy center, analogous 
to the Stark effect in electromagnetism. More precisely, the eccentricity vector e 
of the orbit is the sum of a fixed vector e^? directed towards the Galaxy center 
(proportional to the difference of the accelerations of the two bodies) and of a 
rotating vector eR{t) corresponding to the usual relativistic periastron advance 
at angular velocity w^. Several asymmetric systems of this kind (such as PSRs 
1713-K0747, 2229-f 2643, 1455-3330) happen to have a very small eccentricity. The 
only explanation would be that the rotating vector eR{t) is precisely canceling the 
fixed contribution e^? at the time of our observation: e^^-f ej^(t) « 0. However, this 
is very improbable if the system is old enough, and one can use a statistical argument 
to constrain the space of theories. Moreover, by considering several such systems, 
the probability that they have simultaneously a small eccentricity is the product 
of the already small individual probabilities. The allowed tensor-scalar theories lie 
between the two (approximate) dotted lines. This test is much less constraining than 
the others because the Galaxy is not compact enough to be spontaneous scalarized, 
therefore its scalar charge ac ~ ao is small, and the difference of the accelerations 
oc {aA — Q;s)ac ~ {q-a — ao)ao is small too. 

To detect the gravitational wave signal from an inspiralling binary with the 
LIGO/VIRGO interferometers, one will perform a matched filter analysis of the 
signal, using the gravitational wave templates predicted by general relativity. In 
tensor-scalar theories, the waveforms are very different, because of the extra con- 
tributions of the scalar waves to the energy loss, Eq. I|1HI . Therefore, the ac- 
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tual detection of an inspiralling binary with the general relativistic templates will 
constrain the magnitude of the scalar contributions. The hatched curve labeled 
"LIGO/ VIRGO NS-BH" displays the region of the theory space which would be 
excluded if a lAniQ neutron star-lOm© black hole system is detected with a sig- 
nal to noise ratio S/N = 10 (the excluded region lie on the hatched side). The 
detection of a double neutron star system with masses similar to those of PSR 
B1913+16 would exclude the bubble of theories labeled "LIGO/VIRGO NS-NS". 
As we can see on Fig. 0] these regions are already excluded by binary-pulsar tests. 
Therefore, we can conclude that the general relativistic wave template do suffice 
for these interferometers: Even if there exists a scalar partner to the graviton, we 
anyway already know that it is too weakly coupled to matter to change significantly 
the waveforms. This is a good news, because the inclusion of possible scalar con- 
tributions would have considerably slowered the data analysis. On the other hand, 
it has been prover^l that the LISA space interferometer could still be sensitive to 
scalar-field effects. Indeed, the detection of a neutron star inspiralling a lOOOm© 
black hole would probe values of [--yPP^ — 1| as small as 4 x 10^®, i.e., an order of 
magnitude tighter than the best present bounds Q [or a factor 3 tighter for the 
matter-scalar coupling strength ag]. However, we will see below that the binary 
pulsar PSR J1141— 6545 will probably probe similar values of j^^P^ — 1| around 
2010. Therefore, it may not be necessary to start including scalar corrections to the 
wave templates for LISA: This binary pulsar should tell us, just before the launch 
of the LISA mission, whether there is or not a scalar contribution at this order. 

The constraints imposed by the recently timed PSR J1141— 6545 are much 
tighter than the previous ones, and we plot them as a dot-dashed curve in 
Fig. |S1 We pre sent here preliminary results, which will be refined in a forthcoming 
publication.l^ Note that the vertical scale of this Figure has been expanded by a 
factor 2 with respect to Fig. 01 To ease the comparison of these two Figures, we have 
repeated in Fig.Elthe curves corresponding to PSR B1913-I-16 and to solar-system 
constraints. The reason why PSR J1141— 6545 is so constraining is because this is 
an asymmetrical system, composed of a neutron star and a white dwarf. Therefore, 
like PSR B0655-f64 above, it generically emits a large amount of dipolar scalar 
waves, which are inconsistent with the small measured value of its orbital period 
derivative P. The lower-left panel of Fig. ^ shows that the experimental precision 
on this quantity is not very good yet. However, this is already much better than 
the mere upper bound on P that we knew for PSR B0655-I-64. 

Like all binary pulsars, PSR J1141— 6545 excludes tensor-scalar theories cor- 
responding to Po < —4.5, even for a strictly vanishing ao. This is again due to 
the phenomenon of spontaneous scalarization described in Sec. 3. But it is quite 
remarkable that this system is also very constraining in the region of positive /3o's. 
It is almost as powerful as the previously known solar-system bounds (|2Jl. This 
may seem paradoxical, because we saw that when /Jq > 0, strongly self-gravitating 
bodies are much more weakly coupled to the scalar field than matter in the solar 
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Figure 5. Constraints imposed on the theory space by the two recently timed binary pulsars PSR 
J1141— 6545 and J0737— 3039. The thinner lines display the constraints which will be reached when 
the orbital period derivative P is determined with 1% accuracy. We also display here the impressive 
solar-system constraint obtained recently from the observation of the Cassini spacecraft near 
solar conjunctionl3 



system. The reason why this binary pulsar anyway provides a powerful test in this 
region of the theory plane is again its asymmetry. Indeed, the pulsar's scalar charge 
is exponentially small, \ola\ <SC |q^o|i whereas the weakly self-gravitating white dwarf 
companion has a standard (background) scalar charge as « ao. Therefore, the 
dipolar radiation term (|13ll . proportional to [ua — olbY ~ "^oi iioi^ negligible, 
and the small observed value of P is still constraining. The asymptotic limit of the 
bound on |ao| for /3o +oo can also be estimated analytically,"^ by imposing that 
the dipolar contribution is smaller than the experimental error AP. One finds 



al < 



96 (2'nG{mA+mB) 



5 V 



Pc3 



2/3 



1 + i^^ + i^^ 
(i + f)(i-e2) 



AP 



P 



5 X 10"^ 



AP 



P 



(15) 



With the present uncertainty on P°^^ = (— 4 ± 1) x 10 this gives |ao| < 0.004 



"^Our formal limit /3o — » +oo actually means /3o ~ 10 or 30, but for extremely large values of 
this parameter, the "de-scalarization" phenomenon would affect the white dwarf too, and the 
constraints would thus weaken. 
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for (3q — > +00, consistently with the curve plotted in Fig.|Sl As mentioned is Sec. 2, 
the precision on P should reach 1% by the end of the decade. This would give an 
asymptotic bound |q;o| < 7 x 10~^ for /Jq — > +00, corresponding to |7^^^| < 10~®, 
more than one order of magnitude tighter than the recent limit (jSJ • This estimate is 
confirmed by our numerical plot, in Fig.|3 of the constraints that PSR J1141— 6545 
will impose when this precision is reached. Of course, a detailed analysis of the 
possible sources of noise will be necessary then, notably of the tidal effects on the 
white dwarf, and of the Doppler contribution due to the acceleration of the system 
towards the center of the Galaxy. 

To give the reader a better feeling of the above limits, we display in Fig. 
the mass plane for PSR J1141— 6545 in four different scalar-tensor theories. The 




Figure 6. Mass plane {niA = pulsar, ms = companion) for PSR J1141— 6545 in four different 
tensor-scalar theories of gravity. Only the upper-left plot corresponds to an allowed model (at the 
la level). In the lower-right plot, the system is assumed to be a double-neutron star system. 

upper-left plot corresponds to Brans-Dicke theory, i.e., to a linear matter-scalar 
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coupling function a{ip) = aQip, with Uq — l/(2a;BD + 3)- The value of uq has been 
fine tuned to pass the test exactly at the la level. This theory corresponds thus 
to the intersection of the thick dot-dashed curve of Fig. |31 with the vertical axis. 
The corresponding limit |7^^^| < 6 x 10^"* wbd > 1800) is not far from the 
unpublished VLBI constraint pl as can be also directly seen on Fig.|31 Note that in 
Fig. El the P strip has been significantly displaced with respect to its location in 
general relativity (lower- left panel of Fig. 

The upper-right panel of Fig. displays the mass plane for a theory which was 
allowed by all experimental data at the time of the MGX conference, but which 
is violently ruled out (by 5a) by PSR J1141— 6545. Indeed, the three strips do 
not have any common intersection. This illustrates that this binary pulsar is much 
more constraining than the others, at least in the privileged class of tensor-scalar 
theories. However, the subsequent solar-system bound (0) does even better. 

The lower-left panel of Fig. |5| provides an example of a theory in which the 
three strips 7t, uj and P do have a common intersection, but it is located is the 
shaded region corresponding to | sini| > 1 (see Sec. 2 above). [The lower-right panel 
gives another example.] Therefore, in such a case, this (Keplerian) mathematical 
constraint becomes crucial. If the inequality | sini| < 1 is not taken into account, the 
bump of the thick dot-dashed curve blows up to very large values of |ao|, in Fig. El 
Note however that such theories are anyway ruled out by solar-system experiments. 

The lower-right panel of Fig. illustrates the strange shapes that the curves 
can take in some theories. Here, the topology of the P "strip" has even changed. 
Note also that the intersection of the strips corresponds to masses which are signifi- 
cantly different from those obtained in general relativity (lower-left panel of Fig.^. 
Therefore, it would be inconsistent to use the masses obtained within general rel- 
ativity to compute the predicted P in another theory. Obviously, all the strips 
must be computed within the same theory, and the binary-pulsar test is passed 
if they have a common intersection, even if its location differs from that obtained 
in general relativity. Contrary to the other three mass planes, the lower-right one 
has been computed while assuming that PSR J1141— 6545 is a double-neutron star 
system. This is excluded at the 90% level by the formation scenaricPof this system. 
However, it remains instructive to study the dependence of the constraints on the 
nature of the companion. If it were a neutron star, the system would be much 
more symmetric (a a ~ as), therefore the dipolar contribution (|13|l to the energy 
flux would be much lower, and the observed value of P less constraining. We did 
plot the corresponding bounds in the theory plane (|ao|,/3o)i but we do not display 
them here to clarify the Figures. It suffices to mention than PSR J1141— 6545 would 
give constraints similar to those of PSR B0655-f 64 if it were a double neutron star 
system (see Fig.0J. 

Let us come back to Figure jSl above, where we also plotted as dashed lines the 
bounds imposed by the double pulsar PSR J0737— 3039, at present and when a 
1% accurate P is available. Curiously enough, in spite of the great precision of the 
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measures, and although the mass ratio mA/TiiB = 1-07 is tightly fixed independently 
of the theory, this system appears not to be extremely constraining for scalar-tensor 
theories. When a 1% accurate P is measured, it will of course provide a much 
stronger test, and become more constraining than the Hulse- Taylor binary pulsar 
PSR B1913+16. However, it pales in comparison with the neutron star-white dwarf 
binary PSR J1141— 6545, even with its present large uncertainties on P. The reason 
is that the two components of PSR J0737— 3039 are neutron stars with similar scalar 
charges, and therefore that its dipolar radiation (|13|l is weak. One may also wonder 
why at present, without any observed P, the constraints imposed by this system 
are so loose, as compared to those of PSR B1913-f 16. The reason seems to be that 
the five observed functions of the masses (c/. the lower- left panel of Fig. ^) depend 
only weakly on the presence of a scalar field. Two of them, r and xa/xb^ are even 
totally insensitive to it. Only 7t can vary a lot when spontaneous scalarization 
occurs (because of the blowing contribution as91n/yi/9(^0i see Sec. 3), but this 
post-Keplerian parameter happens to have still large experimental uncertainties. 

Figure [3 displays the mass plane {mA,mB) for PSR J0737— 3039 within two 
tensor-scalar theories. The left panel illustrates that this system has the capabil- 




Figure 7. Mass plane [rriA = pulsar, ms = companion) for PSR J0737— 3039 in two different 
tensor-scalar theories of gravity. The left plot corresponds to a theory which is ruled out by all 
binary-pulsar tests, including the present one. The theory considered in the right panel passes the 
present test, although it is ruled out by several other binary pulsars. 

ity of ruling out some models. As shown in Fig. |S1 it forbids all theories with 
/3o < —4.5, even for a strictly vanishing ao, i.e., even if the theory is strictly indis- 
tinguishable from general relativity in the weak-field conditions of the solar system. 
This property is shared by all binary pulsars. They exhibit similar deformations 
of the various strips in the mass plane, as soon as —(3o is large enough. Notice 
in particular the characteristic shape of the 7t strip, caused by the spontaneous 
scalarization of neutron stars. 
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On the contrary, the right panel of Fig. [7| illustrates that this system is presently 
less constraining than several other binary pulsars. Indeed, the model ao ~ 0.05, 
/3o = -3 is inconsistent with PSRs B1913+16, B0655+64, and J1141-6545, but 
it passes the test for the double pulsar J0737— 3039. This mass plane confirms 
that the five observables are weakly dependent on the scalar field. The largest 
deformation occurs for the 7t strip, with respect to general relativity (lower-right 
panel of Fig. P), but its width is large enough for the test to be passed. Of course, 
such a model will be also ruled out by this double pulsar J0737— 3039 once its P is 
measured with reasonable precision. 

Let us recall that this double pulsar will anyway provide brand new tests of 
relativistic gravity, as well as very important information about the astrophysics 
of pulsars, although its discriminating power seems rather weak at present in the 
framework of tensor-scalar theories. Moreover, if the new tests happen not to be 
passed by general relativity, one can already bet that no tensor-scalar theory will 
be able to pass them either. Therefore, this system has actually the capability of 
ruling out the best class of gravity theories as a whole. On the other hand, if general 
relativity passes the new tests, as we expect, this double pulsar will not teach us 
much about tensor-scalar models. 



5. Conclusions 

Binary pulsars are ideal tools for testing relativistic gravity in the strong field 
regime. In the most natural class of alternatives to general relativity, tensor-scalar 
theories, their probing power has been shown to be qualitatively different from weak- 
field experiments: They have the capability of testing theories which are strictly 
equivalent to general relativity in the solar system. 

Two fantastic binary pulsars have been timed recently. The double pulsar 
J0737— 3039 promises to be the best laboratory for testing general relativity itself, 
and for studying the physics of pulsars. On the other hand, the neutron star-white 
dwarf system PSR J1141— 6545 is by far the most constraining binary pulsar known 
at present, because its asymmetry implies that it generically emits strong dipolar 
gravitational waves in scalar-tensor theories. It is already almost as constraining 
as solar-system tests even in the region of positive /3o's, where binary pulsars never 
competed with weak-field experiments up to now. It should probe values of the 
Eddington parameter [^^^^ — 1| ^ 10^^ by the end of the decade, i.e., more than 
one order of magnitude better than present solar-system limits. 

Binary pulsars are so precise that they already exclude the models which would 
have predicted some significant scalar-field contributions to the gravitational wave 
templates necessary for LIGO and VIRGO. Therefore, one may use securely the 
general relativistic templates for these interferometers. It is still possible that such 
scalar-field effects be detectable with the LISA space interferometer, but binary 
pulsars will probably give us tighter bounds before it is launched. 
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